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ABSTRACT 


A coil antenna consisting of approximately 6000 turns 
of copper wire was utilized to measure the horizontal com- 
ponent of fluctuations of the earth's magnetic field on the 
floor of Monterey Bay in water depth of approximately 50 
meters. The results indicate that the power spectral den- 
sity of the fluctuations varies from 10nT2/Hz SE 
1079һт2/Н2 at 25 Hz, a monotonic decrease of about 6 dB/ 
octave, except in the 8-20 Hz region where the Schumann 
resonances occur. While the sensitivity of the equipment 
was insufficient to measure the vertical component of the 
Mie tuation we can put an upper limit of 1072nT7/Hz at 1 Hz 
and 10-6(nT)2/Hz at 10 Hz on the magnitude of this component. 

The same sensor was also used to measure various com- 
ponents of the field fluctuations at a remote land site 
(Chew's Ridge). In the frequency range observed the gen- 
eral shape of the spectra was similar to those obtained at 
sea. However, a strong azimuthal variation at certain 
frequencies was noted in the land data. The possibility 
that these directional signals are of man made origin 


cannot be excluded at this time. 
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1. INSURODUOTIPON 


This thesis research is a part of an ongoing effort, at 
the Naval Postgraduate School, to obtain improved long term 
data and interpretations of the electromagnetic noise on 
the ocean floor. The overall project objectives are to 
study and interpret signals covering 5 decades of frequency 
meom .001 Hz to 100 Hz through the combined use of total 
field magnetometers and ULF/ELF receivers. The overall pro- 
ject emphasizes the importance of obtaining measurements of 
geomagnetic noise fluctuations on the sea floor and on land 
over a period of several years and at various locations. 

The particular objectives of this thesis research are 
eco collect local magnetic field fluctuation data at a 
relatively quiet land site near Monterey, California for 
comparison with data obtained from the floor of Monterey 
Bay and 2) to test a coil type antenna specifically de- 
signed for deep sea use, as well as the required electronics 


and data recording equipment. 


TE 





IT. BACKGROUND 


A. SOURCES OF GEOMAGNETIC FLUCTUATIONS 

Geomagnetic fluctuations result from sources internal 
and external to the earth's surface. Internal sources, 
which account for approximately one-third of the variations 
in earth's magnetic field, are referred to as secular vari- 
ations, having periods of thousands of years. Geomagnetic 
storms, diurnal variations, and micropulsations are due to 
external sources and have periods on the order of days and 
less. 

Probably the best known sources of geomagnetic fluctua- 
tions are those generated by the interaction of charged 
particles with the magnetic fields of the magnetosphere. 
Figure 1 indicates an approximate power spectrum of vari- 
ations measured on the earth's surface having periods of a 
day or shorter. Although this phenomenon is not fully 
understood, much effort has been expended in studying the 
sources and predicting the occurrences of these variations. 
The high frequency fluctuations, referred to as micropul- 
sations, fall into two general classes: Pi or irregular 
und ол and Pe or continuous pulsations and are be- 
lieved to result from wave- particle interaction in the 
magnetosphere (Jacobs 1970). 

Of particular interest are long and short term predic- 


tions of the Pc 1 geomagnetic pulsations (.2-5 Hz). A 


1 
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number of methods are available for the forecasting of long- 
term (Fraser-Smith, 1972 and Matyeyeva et al., 1972) and 
short-term (Fraser-Smith, 1980) geomagnetic pulsation 
occurances. However, the accuracy of these methods must 
eventually be determined by comparison with the observed 
data for both research and operational reasons. 

Thus a better understanding of the time dependent 
component of the geomagnetic field is of considerable 
interest both in geophysical and naval applications. Pre- 
diction of geomagnetic fluctuations is especially valuable 
in Navy operational and projected systems such as extreme- 
ly low frequency (ELF) communications, submarine detection 
(Magnetic Anomaly Detection), and weapon deployment (mines, 
particle beams, and torpedoes). 

In addition to the information about magnetospheric 
processes which can be deduced from a study of time varia- 
tions of the geomagnetic field, these fluctuationa also 
serve as an important tool in the study of other geophysical 
quantities. Thus both conventional and naval magnetotelluric 
techniques make use of the variable geomagnetic field in 
the measurement of ocean floor conductivities. Similar 
measurements are employed in the study of earthquake faults, 
the location of aquifers and geophysical prospecting in 


general. 
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Bi ELECTROMAGNETIC PROPAGATION IN LAYERED CONDUCTING MEDIA 


To describe electromagnetic propagation in layered con- 


ducting media, a three layer model was selected (see figure 


2). Note that plane wave propagation is assumed. This mod- 


el was used to predict the expected geomagnetic variations 


in the media and at their interfaces. Before proceeding 


with an explanation of the model and its use, the follow- 


ing definitions are presented: 


A x 


The 


Wavelength, meter 

Phase velocity, meter/second 

a+38, Complex propagation constant, vente 
Angular frequency, radian/second 
Permeability, henry/meter 


Dielectric constant, farad/meter, Е== 6 where oa 
107° /36n farad/meter 


Relative dielectric constant, for air e =l and for 


sea water = =81 


Conductivity, siemens/meter, for sea water g 4-4 


siemens/meter and for sea floor 93.1 siemens/meter 
Electric field strength, volt/meter 

Magnetic flux density, tesla 

Skin depth, 1/a, meter 

Frequency, hertz 

Intrinsic impedance, ohm 


media are assumed to be non-magnetic with H]=4,= 


Hien DÉI " henry/meter. Sea water and the ocean floor 
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are assumed to be good conductors in the frequency of in- 
Est (04 to 25 Hz), i.e., D/we>>l (conduction current 
dominates). On the other hand, air is assumed to be a good 
insulator, i.e., G/we<<1l (displacement current dominates). 
For these conditions the following equations for the propa- 
gation constants and intrinsic impedances along with their 
numerical values are valid: 
Region 1: (air) 
=: > » c x 
T4730 (ugs 4) O CIO )f meter 
n,” (HA£ ) In. |=377 ohm 
d E il 
Region 2: (sea vater) 


: 5 -3, kb zii 
To= (joua) ° |1.1=5.6(10 J = 
L 


^ ohm 


n.,= (jan /0,) ° İn, 1-1.4(1077)£ 

Region 3: (sea floor) 

Гәв (Јондаз) * İr, 1-8.9 (1077) £? meter | 

n47 Geug/o 4) ° In, |=8.9 (1073) £” ohm 

For simplicity, only the theory for a normally inci- 
dent, plane polarized wave is considered. The modification 
for the general case of oblique incidence is easily made. 
It can be shown that the direction of propagation in the 
sea is vertically downward, regardless of the angle of 
incidence. 

Use of Maxwell's equations and the above assumptions 


results in the following wave equations for the media of 


interest: 
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Region 1: 
r 


Region 2: 
(Sea Water) 


B- 
d 2. 
Region 3: 
(Sea Floor) 
Et 
t7 
Pleure 2. THREE-LAYER “MODEL 
Region 1: (air) 
E - Ee i!?,g c! 1? 
XI l È 
B E -T1ız Tız 
y1= (ug/n) (Ee *E.e } 
Region 2: (sea water) 
m ei?” 
X2 + - 
= mıza [52 
D (ug/n4) tE,e b e 4 
Region 3: (sea floor) 
= -T32 
ET Е,е 
. -T32 
Ser? (ug/n4) E,e 


The boundary conditions require that the electric and 


magnetic fields be continuous at each.interface, i.e., at 


b 





the air-sea water and sea water-sea floor boundaries. 


These conditions are: 


=f) e + = + 
z=0: E,+E, E, FE_ 


(4y/n)TE-E,) = (u¿/n3)(E,-E_) 
Dee: Se "ap el"? = = ә 
Sg Tod E zit 
(ug/n4) tE,e Ee * } = (ug/ny)E,e 


We have here five unknowns with only four relations 
between them. To obtain a meaningful result from the above 
equations, the various E components are expressed as 


ratios, normalized to Ei. The form of the ratios are: 


E. - (E,/E,), E, = (E,/E¿), E, - (E /E,), and E. = 
(E, /E:) 
Algebraic calculations result in the following expressions: 
2n2 (n2+n 3) 
= =2P24d 
— msena) (nitnz)-(na-n1) (n3-n2)e ^ * 
2n » (nas-n2) 
2852177) (nin )e”” 774 (n2-n:) (n:-nə) 
a -(T2-T:)d (lT2*l'3)d 
E, - Е,е +E e 
BT E +E -1 


Using the above equations for region 2 (sea water), 
an expression for the induced-to-incident magnetic field 


on the sea floor is obtained. 


= SS EZ Бә2 . = 

7: = Е, /na = E, (E,e E e ) since E; (ni/uo)B; 
= "1:0 Гә2 

В, (2) = (П:/ио)В. (Е,е Ee *^) 
m = Rp. I > Z 

B. (z) = Бр. (ni/us) t&, e IN j 


Now assuming that n2 and n3 are much less than ni, we define 


18 





1 
q 7 na/na2 € (02/0317. After some algebraic manipulations, 


we obtain: 


la(d-z) ,. j,,.-Ta(d-z) 
B (z) * 21 (q+1)e – ( 1)е ) 
m? -Г28 


(а41)е129- (а-1)е 
For = and substituting the numerical values of no, ns, 
lUo into the above equation at a specified frequency, one 
can obtain a plot of the induced-to-incident magnetic 
field versus depth. See Figure 3. Clearly, if one is to 
detect magnetic field variations in sea water at any 
appreciable depth, then the variations must be of low 


frequency. 
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C. GEOMAGNETIC FORECASTS 

The Space Environment Services Center located in 
Boulder, Colorado publishes a Preliminary Report and Fore- 
cast of Solar Geophysical Data, including geomagnetic 
activity. A more detailed synopsis of geophysical data is 
contained in the environmental data services monthly publi- 
cation "Solar-Geophysical Data". Of importance to the re- 
search effort here at Naval Postgraduate School are the 
reports of the A and K-indices for the mid-latitude region 
containing Monterey. The major monitoring station for the 
mid-latitude areas in the U.S. is located in Fredericks- 
Викд, Virginia. 

A and K-indices reported by Fredericksburg are of in- 
terest to NPGS, since this station lies approximately on 
the same latitude as Monterey. The A-indices indicate a 
24-hr daily variation in the earth's magnetic field while 
the K-indices indicate a 3-hr period of variation within 
the 24-hr A-index. The K-indices vary in range from O 
(very quiet) to 9 (extremely disturbed) with an A-index of 
30 or greater indicating geomagnetic storm conditions. 

The above information was utilized to correlate geo- 
magnetic predictions and forecasts with geomagnetic data 
obtained in Monterey Bay. This correlation concept is util- 
ized below in Section V, experimental results. (See Appen- 


dix C for a further description of geomagnetic forecasts). 
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D. REVIEW OF EARLIER WORK 

A fair amount of research has been and is being done in 
the detection and prediction of low frequency variations 
(less than 100 Hz) of the earth's magnetic field, both on 
land and in the sea. Earlier experiments by Horton and 
Hoffman (1962), Santirocco and Parker (1963), and Davidson 
(1964) revealed that the power spectral density of back- 
ground geomagnetic activity in the frequency range TO”: 
Hz decreased monotonically vith a slope of approximately -6 
dB/octave. The spectrum then levels off and is dominated 
by the Schumann resonances in the frequency range 1-40 Hz 
(Schumann and Konig, 1954). Measurements indicate a gen- 
eral increase in background activity between 20 and 70 Hz 
as reported by Campbell (1966), with a continued sharp 
decrease in slope above 80 Hz. 

Spectrums obtained by Larsen and Egeland (1968) indi- 
cate that the general decline observed at frequencies below 
5 Hz continues through 70 Hz and probably to higher fre- 
quencies (Maxwell and Stone, 1963; Westerlund, 1970) with 
the Schumann resonance activity being superimposed on the 
declining background as a form of fine structure. Hence, 
there is a remarkably steady decline of the background fluc- 
truations at a rate of approximately 6 dB/octave covering 
almost 6 decades of frequency (1077 - 10? Hz) (Fraser-Smith 
ad Buxton, 1975). 


2° 





Recent measurements conducted by Fraser-Smith and Bux- 
ton in the .1 - 14 Hz frequency range confirmed earlier 
predictions and detections. They discovered that the gen- 
eral decline inb=ekerounda activity continued for frequen- 
cies up to 5 Hz, but for frequencies above 5 Hz the decline 
appeared to be suppressed by the occurence of the Schumann 
resonances. 

More recent work was accomplished in the low frequency 
ranges at La Mesa in Monterey. Barry (1978) and Clayton 
(1979) investigated geomagnetic fluctuations in the fre- 
quency ranges .1-10 Hz and .4-40 Hz, respectively. Their 
results are in general agreement with earlier measurements. 

Chaffee (1979) measured geomagnetic variations of the 
earth's magnetic field on the floor of Monterey Bay in 100- 
200 meters of sea water and covering the frequency range 
.01-3 Hz. Chaffee observed the power spectrum of typical 
geomagnetic fluctuations to vary from 1 nT*/Hz at 0] ELO 
approximately 1075 nT*/Hz. Most low frequency studies made 
in sea water have been performed to detect geomagnetic 
fluctuations produced by surface waves. Since little data 
is available on actual geomagnetic noise at low frequen- 
cies, the project at Naval Postgraduate School is a con- 
tinuing effort to create a reliable data base for future 


reference and correlation. 
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III. EXPERIMENTAL EQUIPMENT AND TESTS 


A. EQUIPMENT CONFIGURATION 
Most of the equipment utilized for this project was 
either 25 and/or manufactured at the Naval Postgra- 
duate School, while test equipment used was mostly of 
commercial origin. Much effort was expended in testing and 
ensuring that the experimental equipment operated as con- 
cieved and predicted. It was important for all equipment 
to be dc coupled to avoid loss of low frequency signals. A 
discussion of the equipment and its configuration follows: 
исә Acquisition Equipment 
SENSOR 
The general hookup of the data acquisition 
equipment is depicted in Figure 4. The sensor is a self- 
supporting, continuously wound, non-center-tapped coil an- 
tenna manufactured from approximately 6000 turns of 18 gauge 
copper magnet wire. The sensitivity of the coil is estab- 
lished in Section III. B. Its weight is approximately 100 
pounds. Its overall dimensions (.15m by .37m) were deter- 
mined by the largest glass sphere commercially manufactured 
by Benthos, Incorporated. Refer to Figure 5. The coil re- 
sistance is 120 ohms and its self-inductance is approxi- 


mately 9.31 henries. 
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PREAMPLIFIER 

Various differential preamplifiers were tested 
and used including the PAR-113 amplifier manufactured by 
Princeton Applied Research Corporation. This preamplifier 
was selected because of its higher signal to noise ratio, 
low power requirements, and compactness. The preamplifier 
itself is packaged as an integral part of the sensor, i.e., 
it is housed in the glass sphere with the sensor. Attempts 
to separate the preamplifier from the sensor resulted in 
additional noise caused by cable motion. 

The front end of the preamplifier utilized a low- 
27 ’stete of the art" differential device (LH 0038 C) 
manufactured by National Semiconductor Corporation. It has 
a single-ended output which is connected to additional 
stages of amplification. The final stage of the preampli- 
fier contains an active filter which is set to a cutoff 
frequency of 20 Hz. The overall preamplifier gain, in- 
cluding filters, is approximately 66 dB. A graph of the 
preamplifier gain characteristics is included as Figure 14. 
A schematic diagram of the preamplifier circuit is contained 
in Appendix A. 

VOLTAGE CONTROLLED OSCILLATOR 

The schematic diagram of the voltage controlled 
oscillator (VCO) is included in Appendix A. The function 
of the VCO is to convert low frequency analog voltages re- 


ceived from the sensor through the preamplifier into a 
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varying frequency, centered at 1500 Hz. The output level of 
the VCO is approximately 0.2 volt (zero to peak) square wave 
of varying frequency which is recorded on an analog cassette 
tape recorder. The conversion factor for the VCO is 1.25 
Hz/mVo - 
REFERENCE OSCILLATOR 
A 2 KHz reference oscillator signal is recorded 
Simultaneously with the VCO signal on a separate channel. 
The schematic diagram of the reference oscillator is also 
included in Appendix A. The primary purpose of the refer- 
ence oscillator is to reduce tape recorder noise (wow and 
flutter) during playback. 
TAPE RECORDER 
The tape recorder employed for this project is a 
portable stereo cassette type (Marantz Superscope, Model 
CD-330), powered by a 6 Vdc battery. The Superscope was 
selected to meet the size constraints imposed by the glass 
spheres. Additional characteristics are: 
(1) Tape Speed: 1-7/8 in/s 
(2) Frequency response: 40 Hz to 12 KHz 
(3) Signal-to-noise ratio: 50 dB 
(4) Wow and flutter: .12% 
POWER SUPPLIES 
Gould Gelyte lead batteries are utilized to meet 
dc power requirements. They are sealed, rechargeable, and 


virtually maintainance free batteries, containing no free 
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electrolyte. Various sizes are used to meet the power 
specifications of all data taking and analysis equipment, 
except test equipment. 
2. Data Analysis Equipment 
' Refer to Figure 6 for the general hookup of the data 
analysis equipment. The tape recorder and power supplies 
(discussed above) were employed in data analysis as well as 
in data taking. A description of the additional analysis 
equipment follows: 
MIXER 
Appendix A contains the schematic for the mixer 
+ 


circuit. The mixer is powered by © 12 Vdc and accepts the 


recorded sensor data and reference oscillator signals at its 





REF. CHANN FREQ. 
“s... TO VARYING 
RECORDER |DATA CHANNE FREQ.| VOLTAGE [VOLTAGE 


CONVERTER 





SPECTRUM 


ANALYZER 





Figure 6 DATA ANALYSIS EQUIPMENT 


inputs from the tape recorder. The difference frequency of 
the mixer is the input to a frequency-to-voltage converter. 


The mixer circuit also contains a low pass filter with a 
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0050077 Trequeney at approximately 1200 Hz. The purpose of 
the low pass filter is to elimate unwanted harmonics of the 
difference frequency. 

FREQUENCY-TO-VOLTAGE CONVERTER (ANADEX PI-375) 

The purpose of the frequency-to-voltage converter 
(FVC) is to convert the difference frequency from the mixer 
into an analog voltage. The varying voltage is an exact 
representation of the geomagnetic fluctuations detected by 
the sensor. The FVC is powered by 24 Vde and has the 
following additional characteristics: 

(a) Input frequency range: 0-2000 Hz 

(b) Output de voltage: 0-10V full scale 

(c) Conversion factor: 4.91 mV/Hz 
The output of the converter is analyzed with a spectrum 
analyzer. 

SPECTRUM ANALYZER (S/A) 

Numerous spectrum analyzers were used for data 
analysis: The Schlumberger Model 1510-05, The Ubiquitos 
Model UA-500, The Mini-Ubiquitos Model 440, and the Hewlett 
Packard Model 3582A. Each of these analyzers has inherent 
advantages and disadvantages revealed during their operation. 

The Schlumberger S/A does not provide for de 
coupling and hence was not adequate for low frequency appli- 
cations. The Ubiquitos has dc coupling and limited low 
frequency capabilities. Also an external display had to be 


provided to obtain a power spectral density read out which 
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would be highly inaccurate. The Hewlett Packard Analyzer 
was capable of de coupling and provided for a digital dis- 
play of the power spectrum. However, it does not have a 
logarithmic scale for frequency. 

The most adequate S/A utilized was the Mini- 
Ubiquitos Model 440. It provides for dc coupling, has a 
digital display of the power spectrum, and has a logarithmic 
as well as a linear frequency scale. 

Since ChewarscriminaLor output is a positive dc 
voltage, it was necessary to provide for de offset in order 
to increase the input sensitivity of the spectrum analyzer. 
The Tertronix Model AM 502 differential amplifier was 
utilized with its gain set to Xl. 

3. Underwater Equipment 
Underwater equipment included the following (see 
Figure 7): 

(a) Benthos vacusealed glass spheres with 
associated hard hats, penetrators, and sealant strip. 

(b) Cable and connectors 

(c) Anchor and nylon line 
The particular spheres employed during this project were a 
low expansion borosilicate type, having a 0.404 meter in- 
Side diameter. The actual spheres used in the tests were 
two hemispherical shells specially matched to guard against 
leakage. They are pressure tested by the manufacturer, to 


1000 atmospheres, which is equivalent to a sea water depth 
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of 10,000 meters. Before the spheres are placed in the 
hardhats, which provide protection against breakage, the 
hemispherical seam is sealed with sealant strip and taped. 

Penetrators were installed in the wall of the 
Benthos Glass Housings to provide electrical connections to 
enclosed equipment. A pressure type penetrator is used to 
activate the tape recorder when the instrumentation sphere 
is at a depth of approximately 10 meters. 

Signals from the sensor to the instrumentation 
sphere are propagated through a twisted-pair, shielded, 
rubber coated, 25 meters long cable with Brantner type sea 
connectors. The connectors are made from a special molded 
neoprene material; they have brass, hard gold plated con- 
tacts and are pressure rated to greater than 1000 atmos- 
pheres when mated. 

An anchor supports the sensor sphere and provides 
added negative buoyancy to prevent the system from floating. 
It was designed and constructed from PVC piping, filled 
with sand, and neoprene-coated lead bricks were strapped to 
the piping for additional weight. Extreme care was taken 
to ensure that any exposed metallic surface which might 


induce unwanted currents into the sea water are eliminated. 
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PSR GENSOR SENSITIVITY 
1. Theoretical Sensor Sensitivity 
The use of wire coils to measure magnetic field 
changes is not new. Faraday's Law of Induction states that 
the electromotive force (emf) induced in a single-turn 


circuit is equal to minus the rate of change of flux: 


emf =- (3-1) 


Viv 
cr ho 


The flux is that which passes through a surface 


bounded by the conductor and is defined as: 
6 = fB-da (3-2) 


If we consider N turns of wire connected in series, 
then the total induced emf is the sum of the emf induced in 


each separate turn, or for a loop: 


90 


Jr us) 


emf - -N 


In this work it shall be assumed that magnetic 


fields vary sinusoidally with time and may be written as: 


5 = Be Jut 
O 


(3-4) 
Where B is a magnetic vector field and only a 
function of spatial coordinates and not of time. By substi- 


tuting equation (3-2) into (3-3) and with equation (3-4) we 


have: 


emf = -NB JVA (3-5) 
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where A is the area enclosed by the coil. The magnitude of 


this emf is; 


eme NB. (27f)A (3-6) 


If the.output of the sensor is connected to a high input 
impedance amplifier, then negligible current can flow 
through the sensor which might otherwise produce a back emf 
and reduce the magnitude of the field to be measured. The 
sensor's open circuit voltage or sensitivity can be calcu- 
lated from the parameters given in Section III. A. 1 fora 
5007 eld at a frequency of 1 Hz; emf - 3.1117 uV, thus 20 
log (emf) = -110.14 dB. (re: 1 Volt) 

The area used in calculating the emf for a single 
turn of wire is easier to visualize than the area of a 
finite width coil. It was decided that an average area 
obtained from the area within the inner radius and the area 
within the outer radius would be used. This method appears 
to be more accurate than using an average radius to compute 
the area. Figure 8 depicts the expected sensitivity of 


the sensor. 
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2. Experimental Determination Of Sensor Sensitivity 


Two methods were utilized to determine the sen- 
sor's actual sensitivity. The first method consisted of 
producing a known field at a specified distance from an 84 
turn copper coil and measuring the electromotive force or 
voltage signal produced by the sensor. The second method 
was similar to the first, but a .61 meter radius Helmholtz 
coil was used to establish a more uniform magnetic field in 
the central region. 

METHOD I. 

The magnetic field produced by a circular ring of N 

turns at a point Z along the symmetry axis is given by: 


2TNIr? (10?) 3 


LE) = = b. 57 
e + yr?) /2 


Equation (3-7) is obtained from the Biot-Savart Law using 
the units of I (ampere), r (meters), 2 (meters), and B 
(nanotesla) as depicted in Figure 9. Test coil specifica- 
tions were as follows: 


O 
84 turns 


E 


N 


r inner = 15 cm 


r outer 16.5 cm 


where r was computed from the average area of the coil. 


Rearranging equation (3-7) we have: 


3 
(22 + x?) 7 B. (3-8) 


eoo rn 


Irms = 
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Specifying the desired field strength and distance deter- 
mines the current through the coil at different frequencies. 
To produce the required current a 1102 resistor was placed 
in series with the coil of 42 depicted in Figure 10, and 
the Basen voltage was set according to: 


s = . + R). (3-9) 


—Л 


The current was determined to be within * 1.04€ of one mA and 
the field produced along the central axis of symmetry was 
determined to within 1 10.% of a nf. 

The data obtained on the sensor's sensitivity for a 
5 nT and 66.4 nT signal, when compared to the theoretical 
curve, revealed that a more uniform field in the vicinity of 
the sensor was needed and that the test should be performed 
in an area remote from metallic structures such as ina 
laboratory. The data obtained was within 6 dB of the theo- 
retical values as can be seen in Figure 12. 

METHOD II. 

The second method, that of using a pair of Helm- 
holtz coils, produced a more uniform field over a larger 
area. Figure 11 illustrates how the coil was designed to 
produce an easily accessible, almost homogeneous magnetic 
field in the central region between the coils. 

Expressions for the magnetic induction vector for 


the axial and radial componets are: 


oY 








Iı mü 3 zr” 

000 er 25-.£ç2 SE 3-10 
MA TR AA 
pee r SDI HE) 27577.) (3-11) 
m S/5R 

as given in (Jefemenko, 1966). 
Senin the central region: | 
B = He (3-12) 
By oR 
B. = 0 5517 


to the terms of the order (5207 where L is the linear di- 
mensions of the region. Using the same method to determine 
the current in the loops, as explained above, the data 
shown in Figure 13 indicates that the actual sensitivity 
was about one dB less than the theoretical value for the 
entire frequency band of interest. It was determined that 
the accuracy of the measured sensitivity was well within 
the expected theoretical limits. The major source of error 
in the theoretical computuation was determining accurately 
the area enclosed by the conductor due to the finite width 


She coil. 
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Figure 10 Hookup for Establishing a Known Field 





Figure 11 Helmholtz Coil Schematic 
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To obtain the corrected power spectral density referenced 


ipo 1 nT^/Hz, the following procedure (derived in Appendix 


D) was used: 


N dB 
(A) 


(B) 
(C) 


(C) 


The Spectrum Analyzer reading in dB volts. 


The Hanning window bandwidth correction for 
the particular Spectrum Analyzer used. 


ипи ers iimersom deraned above. 


The corrected power spectral density where 
O dB is referenced to 1nT*/Hz. 


N dB - (B) - (A). 
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C. SYSTEM TRANSFER FUNCTION 

The system can be represented by a simple block diagram 
that converts a varying magnetic field into a voltage that 
can be displayed and whose power spectral density can be 


measured with a Spectrum Analyzer: 


ikən SYSTEM X number of volts 


The transfer function is frequency dependent. It was cal- 


culated by separating the system into its components as 


follows: 
(v) (V) (Hz) (Hz) 
1 nT SENSOR PREAMP | VCO AMP 
X VOLTS 
PVC —— a 


— — — —. 


Once the sensor sensitivity was ти using the Helm- 
holtz coils it remained to compute the gains and conver- 
sions of the electronic components to determine the overall 
system transfer function. 

The preamplifier was designed to be a low pass, low 
noise amplifier with a very sharp cutoff frequency at 20 Hz. 


This cutoff frequency is necessary to attenuate the 60 Hz 


7” 





fields present at all land sites. As can be observed in 
Figure 14 the gain characteristics are flat to 10 Hz with 

ии Gf 600.5 dB. The sharp cutoff at 20 Hz was obtained 
With a final stage active filter network which added slight- 
m 5 the gain prior to cutoff. 

The voltage controlled oscillator (VCO) was designed 
to have a center frequency of 1500 Hz with an output voltage 
level of .4 volts (Peak to Peak). The purpose of having a 
low output voltage at 1500 Hz is to minimize crosstalk be- 
tween tape recorder channels and to operate well within the 
frequency response of the tape recorder. Since the tape 
recorder and mixer (described in Section III. A) have unity 
gain, the overall gain of the VCO and frequency-to-voltage 
converter (FVC) was measured by amplifying the 1500 Hz 
signal to the required level to drive the FVC with a low 
noise, PAR-113 amplifier. The low pass gain characteristics 
of the VCO/FVC combination can be seen in Figure 15. The 
response is similar to that of the preamplifier and is 
characterized by a constant gain of 34.7 dB from de to 4 Hz 
and a slight decrease after 10 Hz as required for additional 
ee S filtering. 

The system transfer function was calculated by adding 
the gain of the preamplifier to the gain of the VCO/FVC com- 
bination and to the sensor sensitivity at equivalent fre- 
quencies. This can be seen in Figure 16. The data presented 


in this research is referenced to 1 nT rather than 1 volt. 
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D. SYSTEM NOISE 

The best amplifier performance is realized when the 
amplifier causes the least decrease in the overall signal- 
to-noise ratio. In many instances the thermal noise gener- 
ated by the signal source resistance is the dominant factor 
in determining the input signal-to-noise ratio. Amplifier 
noise performance can be specified by the amount of noise 
the amplifier adds to the amplified source thermal noise. 
Expressed in decibels, and defined as the "noise figure", 


Gere NF, it is given by: 


NF = 20 log total rms noise voltage at amp output (3-14) 
gain x (source thermal noise voltage 


where the Source Thermal Noise = y4KTRaf volts rms, and 


|I«NBolkbemann's constant, 1.39 x 10742 joules/9K 


D absolute temperature in degree Kelvin 


f = Equivalent Noise Bandwidth in Hz 
R = source resistance in ohms 
The total output noise may be converted to an equivalent 


input noise by dividing by the amplifier gain. Thus the 


noise figure becomes: 


NF = 20 log total rms noise voltage referred to amp input 
source thermal noise voltage (rms 
дә) 
The equivalent noise bandwidth (Af) used in determining the 


source thermal noise is simply 7/2 times the high fre- 


gencv rolloff. 
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The system's preamplifier has a 3 dB point on the low 
pass gain curve at 25 Hz. This yields an equivalent noise 


bandwidth of 
250005 7 - 39.27 Hz 


and a source thermal noise of 8.69 x 1077 rms volts, vhich 
is -161 dB re: 1 volt. 

A 1200 resistor vas connected across to the preamplifier 
input vhile the rms voltage vas measured on the spectrum 
analyzer, By dividing the total output noise by the amp- 
lifier gain the Noise Figure at one hertz was obtained as 
15 dB (see in Figure 17). 

Once the total rms noise voltage referred to the pre- 
amplifier's input was determined, the noise introduced by 
the remainder of the system was determined to be negligible 
in comparison to the preamplifier noise, with the exception 
of four very discrete frequency components introduced by 
the tape recorder. These are centered at 3.4, 6.1, 6,8, 
and 8.1 Hz and have been deleted from the data curves. The 
purpose of the reference signal and mixer is to eliminate 
the tape recorder noise. However, the noise reduction 
is limited to about 10 dB at these four undesired frequen- 
cies. 

The total system noise rms voltage level was measured 
several times at various locations. In particular, it was 


measured each time data was taken in the environment pre- 
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sent (i.e. land or sea). On each occasion the total noise 
was found to be at the level shown in Figure 17 which was 
exclusively front end preamplifier noise. No other sources 


ub roduced noise above this level. 
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IV. EXPERIMENTAL RESULTS 


A. UNDERWATER DATA 

Of the six hours of total recorded data, half was ob- 
tained at sea and half at a remote land site near Chew's 
Ridge in California. The underwater data was taken at two 
stations designated as Station #1 (36°-38'N, 121°-52'W) and 
Station #2 (36°-39'N, 1219-53"N). The depth, local time, 
date, and frequency band of the measurements are annotated 
on each figure. The sea floor at both stations has less 
than 5° of slope and consists mainly of smooth sand. The 
sensor was arranged so that the central axis of symmetry 
was parallel to the sea floor at an arbitrary angle refer- 
enced to a North/South direction. 

eure 18-19 represent one half hour of recorded under- 
water data taken in a moderately heavy sea state 10-15 
swells) on 14 April 80. The remainder of data taken on 14 
April was not usable due to loss of electrical continuity 
between upper and lower spheres. Figures 20-29 represent 
2.5 hours of recorded underwater data taken in calm seas 
Como April 80. 
B. LAND DATA 

All three hours of land data was recorded on 1 May 80. 
The location was chosen because of its remoteness and lack 
of external power sources or power lines. The position of 


the sensor remained fixed throughout each data run, and 


De 





the orientation was as follows: 


VERTICAL 


SENSOR 


NORTH 


I = 80° at 


3” (Monterey, Calif) 


Figure £ 9 2 
201 o. 180° 
32 0 NA 
33 90° 90° 
34 90° 0° 
35 eo S 
36 309 1809 
37 909 90° 


The only prominent signals observed at less than i Hz were 
on the first run. All subsequent runs revealed preamplifier 


U01Se€ to be dominant at less than 1 Hz. 
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6% DISCUSSION AND CORRELATION OF DATA 

Many variables influence the level of geomagnetic acti- 
vity on the sea floor, thus, an unambiguous comparison with 
the sparce previous data is almost impossible. Three at- 
tempts at sea of three hours duration each were made with 
the coil's axis of symmetry directed upward in order to mea- 
sure the vertical field component. The data gathered indi- 
cated that the geomagnetic activity in the vertical plane 
was below the noise level of the system. This agrees with 
theory which predicts that incident plane waves are refracted 
normal to the surface leaving only a component of the field 
in the horizontal plane in the sea. The data indicates that 
a small verticle component in the sea may still be present 
but was observed to be below preamplifier noise. Thus an 
upper limit on the vertical component of magnetic field 
fluctuations can be established as 1072nT“/Hz dim ar ıq 
107Önn2/Hz at 10 Hz. 

Deala difficult to compare the fluctuations of the 
total magnetic field as obtained from a total field measuring 
system and a system which measures only one component of the 
field. Earlier data obtained with a Cesium magnetometer 
(Chaffee 1979) on the floor of Monterey Bay indicated that 
the 1 Hz fluctuations as measured ina 5 hour period were 
-30 dB referenced to 1 nT2/Hz at a depth of 70 meters. The 
average value obtained with the coil on 18 April at 1 Hz was 
-28 dB referenced to 1 nTZ/Hz at a depth of 45 meters. This 


appears to be in very good agreement. 
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The most interesting frequency range observed was above 
1 Hz. The level of geomagnetic activity was observed to be 
quiet between 1-6 Hz but became active at 9.6, 11.9, 16.5, 
and 19 Hz for the sea measurements. The land measurements 
undicatedusimilar characteristics with predominant fre- 
@iemeies centered at 7.3, 10.7, and 24.9 Hz. Polk and 
Fitchen (1962) made measurements of the N-S, E-W horizontal 
magnetic components and showed that the fundamental mode of 
Schumann resonance is not constant in frequency with time. 
They were able to determine the spectral frequency range as 
7 Hz to 11 Hz. The spectral frequencies at 7.3 Hz and 9.6 
Hz as measured on land and in the sea respectively, corre- 
late well with the fundamental Schumann resonance frequency. 

The spectral frequencies of 10.7 Hz and 11.9 Hz as 
measured on land and in the sea respectively have the typical 
time development for Schumann resonances as reported by 
Polk and Fitchen. The second resonant mode should occur at 
ilk Hz as observed by Balser and Wagner (1960). Since the 
data in this research does not correlate well in frequency 
as a second resonant mode, but does indicate typical wave- 
shapes for Schumann resonance, furthur investigation is 
recommended. 

The 7.3 and 10.7 Hz spectral lines appeared to have 
different polarizations as different coil orientations pro- 
duced varying levels of either frequency. This observation 


could not be established on the sea floor, as no control 
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over the angle with respect to North/South was available, 
but was quite evident in the land measurements. The 24.9 

Hz spectral line was observed to increase toward the after- 
noon and was apparently independent of orientation. Fur- 
thur investigation of these frequencies is required to de- 
termine polarization as the methods utilized for determining 
orientation were inaccurate (* 5%) put the data was still 
reproducable and quite convincing. 

The sea floor data taken on 14 April was higher in 
magnitude than the data taken on the 18 April. In part this 
can be attributed to a higher A-index (9) on 14 April com- 
pared to an A-index (4) on 18 April. Internal and surface 
wave motion was higher on 14 April as evidenced by separa- 
tion of the cable. No wave height correlation could be 
made, since a wave height measuring device was not avail- 


able in the vicinity of the data gathering region. 
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V. EQUIPMENT/SYSTEM IMPROVEMENTS AND RECOMMENDATIONS 


——— —— — Rr M aáÁ— 0 EM D E RD aquqa amam: —w—w—— 


The system utilized for this project was well suited 
to a shallow sea water environment (200 meters). The 
system and equipment have been demonstrated to operate as 
conceived. To deploy this system in a deep ocean environ- 
ment will require several modifications. 

The following discussion on system improvements is 
presented as recommendations to make the system more reli- 
able and usable in the ocean environment. The tape recorder 
is probably the single, major source of equipment troubles 
in the system. Although the record heads and capstan drive 
were meticulously cleaned and demagnetized, the recorded 
information still appeared noisy and distorted. The mixer 
circuit was able to attenuate and filter some of these 
distortions, but not to the degree that totally eliminated 
them. Initial investigation revealed that the longer the 
tapes, the more prominent the noise. Shorter tapes were 
utilized with some success. Further investigation revealed 
that the single capstan drive motor and associated circuitry 
were of inferior design and failed to adequately compensate 
for the changing tape reel speeds. As discussed in a pre- 
vious section, the distortions were finally determined to 
be present at certain discrete frequencies. 

Another obvious requirement for a more reliable re- 


corder is the reduction of recording time brought about 
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£7 7. Snorter tapes. Even though digitizing the data 
would increase data gathering time, the present recording 
system's operation in this role is highly questionable. 
Therefore, it is recommended that this tape recorder be re- 
placed by a more reliable instrument. 

Coils used as sensors require less power to operate 
than the more complicated total field measuring magneto- 
meters. However, coil type sensors have lower sensitivities 
(sensitivity is proportional to frequency) and are unable 
to measure the dc component of the earth's main field. The 
present coil size is restricted by the largest glass sphere, 
as discussed in Section III. The sensitivity of the sensor 
can be increased by increasing the area and/or the number 
7 00 of the coil. However, this Ts not possible in the 
present system. Thus the sensor is limited in its capa- 
bilities to detect geomagnetic fluctuations because of its 
lack of sensitivity at low frequencies. 

In addition, higher frequencies are attenuated more 
than lower frequencies in sea water (see Section II). 
Hence, the use of a coil in a deep ocean environment is 
questionable but not completely ruled out. Modifications 
to the preamplifier (to be discussed below) would increase 
the sensitivity of the system and, possibly make it quite 
feasible for deep ocean use. However, it is recommended 
that continued use of the present system be restricted to 


Shallow water deployment or the system be redesigned for 
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land measurements. That is, a larger sensor be designed and 
constructed for increased sensitivity and used to gather 
geomagnetic fluctuations on land. If three such sensors 
were produced, then all components of the earth's magnetic 
field EC could be measured. In addition, the 
Bolarıization of the horizontal component could be further 
investigated by connecting two such coils and determining 
the coherence between the two measurements. 

State of the art devices limit the design of a pre- 
amplifier which operates at low frequencies. In the present 


system, the noise figure of the preamplifier is 15 dB at one 
hertz. More recently an improved preamplifier has been 


designed and built by Dr. Alan Phillips whose noise figure 
is about 10 dB better. This improvement will greatly 
enhance the sensitivity of the system and make its use in 
a deep ocean environment more feasible. 

In the present system, no means exist for determining 
the horizontal orientation of the sensor axis when deployed 
in the ocean. Observed data indicates that the horizontal 
component of geomagnetic field variations may be polarized. 
(See Section IV on discussion of data). It is recommended 
that an aiming device be designed and constructed so the 
exact orientation of the sensor is known when deployed in 
the ocean. 

An improvement in the splicing of cable and connectors 


is recommended. Presently the Research Department at 
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the Naval Postgraduate School splices the connector pigtails 
to the cable and sphere penetrators. This method was 
satisfactory for shallow water employment, however, a more 
E able method is required for deep ocean usage. It is 
recommended that cables, connector pigtails, and penetrators 
be commercially spliced for greater reliability. 

Presently the only means of data correlation are with 
past data collected at the Naval Postgraduate School and 
the indices (see Section II) reported by Fredericksburg, 
Virginia. Since surface and internal ocean waves are known 
to be a source of field fluctuations, it is recommended that 
a wave gauge (various types are available) be purchased and 
utilized as another means of correlating geomagnetic 
variations. 

It is recommended that the various electronic component 
chassis's be integrated into one or two units and a more 
practical method of connecting power be provided to en- 


hance the system's compactness and increase its reliability. 
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APPENDIX A 
EQUIPMENT SCHEMATICS 


Numerous preamplifiers were tested and utilized during 
the extent of this project. Figure 38 was the most recent 
"State of the Art" amplifier which provided the quietest 
operation, i.e., greatest signal-to-noise ratio. 

Figure 39 is the schematic for the mixer circuit used 
at the output of the analog cassette recorder for data 
analysis. Figure 40 is the 2 KHz oscillator reference 
signal generator used in the analog recording system to 
reduce tape "wow and flutter". Figure 41 is the VCO unit 


used in this experiment. 
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APPENDIX B 


EQUIPMENT USAGE 


Equipment utilized to measure, record, and analyze 
magnetic field fluctuations requires a certain amount of 
maintenance and handling. Below is a suggested list of 
hints which might help prevent equipment damage and pro- 
long its usefulness. 

A. Benthos glass spheres - much care should be exer- 
cised when handling the hemispherical shells because they 
are highly susceptible to chipping around the ground glass 
edges when exposed. The shell halves are serialized and 
appropriately marked for proper mating. When exposed to a 
deep water environment, the spheres are subjected to lower 
temperatures. Retreaval from the ocean depths results in 
the condensation of water. A heat gun may be necessary to 
dry the surfaces and return the spheres to ambient temper- 
ature. Temperature differences may also cause the glass 
housings to become vacuum locked. Evenly applied heat 
about the sphere surface will assist in releasing the 
vacuum. 

The strip sealant used for sealing the shells 
should be applied on a meticulously cleaned surface at or 
near room temperature. Failure to properly seal the 
spheres will result in leakage and possibly, equipment 


damage. Toulene and alcohol are the recommended glass 
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55 nme saluwıons. Caution: Toulene is both toxic and 
poisonous. It should be applied with rubber gloves in a 
well ventilated area. 

B, Electronic Equipment - of particular concern here 
is proper application of power supply voltages. Improperly 
applied voltages may result in integrated circuit and/or 
printed circuit board damage. The suggested application of 
dc power is to first, connect the non-common lead (s), 
being careful not to short the common lead against the 
equipment chassis, and finally, connect the common lead. 
Circuit power should not be applied unless there is a load 
on the input. 

Co Portable Cassette Tape Recorder - The tape recorder 
heads and capstan drive should be cleaned and demagnetized 
after each recording. Failure to do this may result in 
noise spikes that could ruin a hard days work. It is also 
recommended that all tapes be "fast forwarded and rewound" 
before being utilized for data taking. High quality 60 
minute tapes are recommended, since longer tapes appear to 
impose greater drag on this type recorder. Before sealing 
the instrumentation sphere, ensure the record button is 
depressed and the record levels are properly adjusted. 

D. Voltage controlled oscillator (VCO) - the center 
frequency of the VCO should be adjusted to 1500 Hz prior to 
sealing the instrumentation sphere. This can be accomplished 
by adjusting the center frequency poteniometer with the input 
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E. Mixer - an amplifier stage is provided for the 
sensor and reference oscillator (clock) inputs to the 
mixer. With inputs from the clock and sensor, adjust the 
gain poteniometers for a 6-10 V peak-to-peak output dif- 
ference frequency. A 6-10 V peak-to-peak output is 
necessary to drive the frequency-to-voltage converter. 
Verify that the output difference frequency is what it 
should be. Note that this procedure may have to be per- 
formed prior to each data analysis, since record levels 
may vary from one recording to the next. 

F. Gelyte batteries - the batteries should be charged 
prior to taking data. Each battery is labeled and a record 
of charge is maintained in the laboratory. Presently there 
are 2 ampere-hour ratings being utilized in the system. 

The recommended charge procedure is: 

18 ampere-hour: charge at 1 ampere for approxi- 
mately 20 hours for full charge. Do not overcharge for more 
than 24 hours in excess of the ampere-hour rating. 

2.6 ampere-hour: charge at 0.2 amp for approxi- 
mately 15 hours for full charge. Do not overcharge for more 
than 24 hours in excess of the ampere-hour rating. 

Go Connectors/cable: before and after each use in a 
corrosive environment, all connectors should be scrupulously 
cleaned to prevent deterioration. Alcohol and cotton-tipped 
applicators should be used. Further, all male connectors 
should be lightly lubricated with silicon grease to ensure 
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an adequate seal and prevent water intrusion. The locking 
sleeves on the connectors are there to prevent stress and 
are not part of the pressure system. The connectors are 
adequately marked and care should be exercised in mating 
them. The cable and connector pigtails should not be sub- 
jected to undue strains and sharp bends. Intermittent 
and/or complete signal loss may result as well as cable 


separation. 
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APPENDIX C 


Geomagnetic Indices and Forecasts can be a valuable aid 
in data correlation and interpretation. In this Appendix 
we include an explanation of the type of geomagnetic in- 
tion which is available on a weekly basis from NOOA 
(Boulder, Colorado) and which we have found to be very 


pUl in our research. 
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DESCRIPTIVE ENT 
Contents of 


Prelimitary Report" and Forecast of 
Solar-Geophysical Activity 


1 May 1980 


The Preliminary Report and Forecast of Solar-Geophysical Activity 
(PRF) is composed each Tuesday and issued the next day by the joint 
National Oceanic and Atmospheric Administration, Space Environment 
Laboratory (NOAA)/Air Force Air Weather Service (AWS), Space 
Environment Services Center (herein denoted as SESC), Boulder, 
Colorado 80303. It is a continuation of the series of weekly 
reports which began in 1951 and were issued jointly by the High 
Altitude Observatory and NOAA and its predecessors. The current 
series began with serial number PRF 001 on 09 September 1975. 


This is a preliminary report based on data available at the time 
EXESUubIjcation and, tnerzerore, cannct reliaoly be cite or reference 
purposes. It is intended ror rapi sfrisuc..ongsnu or primary use 
to commercial and real time operations. Many of the data contained 
here are subject to later revision or refinement. There are no 

Puy rent or other restrictions in photocopying this report. Archival 
quality data, suitable for more extensive studies, are published in 
the monthly publication "Solar Geophysical Data" (SGD) available 

from the National Geophysical and Solar-Terrestrial Data Center 

777:0 Data Center A), Environmental Data and Information Service, 

D6 NOAA, Boulder, CO 80303. At a still later date, definitive 
publication for much of the data is made in the IAU Quarterly 

Bulletin of Solar Activity. Special reports are also issued from 

time to time by the World Data Center for Solar-Terrestrial Physics. 
















This report is made possible through the combined efforts of 
several groups providing data on a fast time basis. A complete 
list of data sources is given in Note 4 at the end of this descrip- 
tive text. Comments concerning the content of or possible improve- 
ments in the Report are invited. 


Frontisviece 


The frontispiece contains the title of the publication, publication 
number and date, and a photograph symbolizing man's dependence upon 
the radiative output of the sun. The photograph is a frame from 
the time-lapse 35 mm patrol movie taken with the hydrogen-alpha 
patrol telescope of the Space Environment Services Center in Boulder, 
Colorado. The photograph was taken December 9, 1968 by Joe Sutorik, 
chief solar observer, in collaboration with student observer and 
mountaineer Galen McFadyen. The two observers carefully charted 
the positions of the setting sun along the profile of the rugged 
Flatiron Mountains three miles west of the observatory and antic- 
ipated the date when the sun would set behind an accessible rock 
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outcropping. On the selected date, Mr. McFadyen scaled the mountain 
and positioned himself so that his full figure would appear silhou- 
etted against the sun. 


First page 


The first page summarizes the highlights of solar-geophysical 
activity during the past seven days and contains forecasts o- 
expected activity during the next 27-day interval. 


in the activity summary and forecasts, solar active regions are 
identified by the SESC-assigned region number, latitude, Carrington 
(heliographic) longitude, and modified Zurich region classification/ 
sunspot area in millionths of the solar hemisphere on the date of 
maximum sunspot area. 


Significant solar activity is discussed in terms of the character- 
Ҹели of the region of origin, x-ray classification of flares (C, 

17 or X), optical importance o£ flares (SN, 1B, etc.), radio emission, 
energetic particle emission, and geophysical effects. The character- 
istics of near-Earth proton events as detected by both satellites 

and ground-based sensors are discussed. Geomagnetic storms and 
disturbances are described. Observations of visual aurora are 
described when available. Whenever feasible, solar and zeomagnetic 
activity will be summarized and forecast within the framework of 

the following routinely used terms: 


Salar activity: 


Very low - usually only quiet regions on the solar disk 
and no more than five of these; fewer than 
ten C-class subflares without important 
centimetric radio bursts or sudden ionospheric 
disturbances during a 24-hour period. 


Low - usually more than five but less than ten 
quiet regions on the solar disk; only C-class 
subflares without important centimetric radio 
bursts or sudden ionospheric disturbances 
during a 24-hour period. 


Moderate - eruptive regions on the solar disk; fewer 
than five M-class x-ray events with important 
centimetric radio bursts or sudden ionospheric 
disturbances during a 24-hour period. 


High - active regions on the solar disk; several M- 
class x-ray events with important centimetric 
radio bursts and strong sudden ionospheric 
disturbances and/or one or two importance 2 
chromospheric flares or X-class x-ray events 
during a 24-hour period. 
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Very high - one or more regions on the solar disk capable 


of producing proton events; 


one or more 


chromospheric flares of importance 2 or 
greater with outstanding centimetric radio 
bursts (500 flux units above background, or 
and major 
sudden ionospheric disturbances observed or 
expected within a 24-hour period. 


greater), 


X-class x-ray events, 


The x-ray classification of solar flares was initiated by the Space 
Environment Services Center on Ol January 1969. 
attempts to classify solar activity by its energy output, and 

offers at least two distinct advantages compared to the standard 


optical classifications. 
ical significance Of a solar event, 


anis eclassitteacion 


It gives a better measure of the geophys- 
and it provides an objective 


means of classifying almost all geophysically significant activity 
regardless of its location on the solar 5 or az (or beyond) the 


Solar limb. The definitions of class C 


and X events are based 


upon the peak x-ray flux in the 1 to 8 er wavelength range, 


as follows: 


Class Peak 1-8Ä flux 
(o in watt m 4) 
C 57 
M 100 2 77- 
X 255 


For descriptive purposes, 
appended to these letter designations. 
multiplier. For example, 


(6 in erg BUR sec 
-2 


ф = 10 


1072 


¢ > 10 


peak flux in the neighborhood of 3x10 6 watt m 
event would indicate an x-ray burst with peak flux near 5x10 ^ watt 
m ^. This expanded classification is generally used only for 

forecasts are usually issued in terms of the 


description of event 


S; 


broad C, M, and X categories. 
events smaller than 1x10 6 watt m 


Geomagnetic activity: 


-1, 


x € 1071 
-1 


a one digit number from 1 to 9 may be 
The letter then acts as a 
a C3 burst indleacess 


n x-ray burst with 
a an X5 


CO designation may be used for 


Im the definitions of the following adjectives used to describe 

A refers to the 24-hour A-index observed at a 
mid-latitude observatory such as Fredericksburg, VA; not the planetary 
A-index, Ap. The K-indices are, likewise, mid-latitude values. 


geomagnetic activity 


Quiet 
Unsettled 
Active 
Minor storm 


Major storm 


— 
— 


N 


A <7, 


15 xc 157 


19755 
30 <A 
557: 


usually no K-indices > 2. 


usually no K-indices > 3. 


< 30, a few K-indices of 4. 


< 50, K-indices mostly 4 and 5. 


some K-indices 6, 
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or greater. 








Sudden commencements are indicated by beginning times given to the 
nearest minute. Gradual commencements are indicated by beginning 
tines to the nearest hour. 


Immediately following the summary of activity and the forecasts is 
a computer-generated display of classical solar-geophysical indices 
for the dates indicated. The column headings are defined as 
follows: 


OPie2600 MHz: The Ottawa NRC 10 centimeter (2800МН2) radio flux 
observed at approximately 1700Z on that date. 


poP NBR: The official SESC sunspot number reported for that date in 
the Joint USAF/NOAA Solar Region and Activity Summary issued at 
02002 the next day by SESC. Sunspot reports are received daily 
from as many as half a dozen observatories. The SESC staff 
reviews these reports and forms a composite picture of each 
individual region, taking into account such factors as the time 
of observation and the quality of seeing. The composite informa- 
tion for all regions is then incorporated in the 0200Z message. 
The sunspot number is computed according to the formula R-10G-*N 
where G is the number of sunspot groups (regions) and N is the 
total number of individual spots comprising all groups. A 
sunspot number of zero indicates there were no visible sunspots 
on that date, while a blank indicates that no observations 
were possible. 


Fredericksburg and Anchorage A- and K-indices: The daily 24-hour A- 
index and 3-hourly K-indices from the Fredericksburg VA (mid- 
latitude) and Anchorage AK (high-latitude) stations monitoring 
the Earth's magnetic field. K-indices range from О (very 
quiet) to 9 (extremely disturbed). A-indices range from 0 
(very quiet) to 400 (extremely disturbed) with an A-index of 
30 or greater, indicating geomagnetic storm conditions. 


Additional pages 


Additional pages are included in the Preliminary Report as solar 
activity requires, and are comprised of computer-generated displays 
which are described below in their normal sequence. The sequence 
and page numbering may vary occasionally depending upon the level 
of activity, reports of special events, or the inclusion of data 
for past days which could not be published at the usual time. 


Alerts Issued This Week 


A summary of all events which exceeded the SESC alert thresholds 
during the previous week. Listed are the dates, times and the 
activity that resulted in the issuance of the alert. This summary 
will enable alert customers to have a complete listing of alerts 
and may adequately fill the needs of some customers currently 
contacted by telephone. 
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Weekly Energetic Event Summary 


72 : xy 3011) energetic (x-ray intensity greacer than or equal 
to 1x10 2 watt m ^) events occurring during the seven days preced- 
ing the Report composition date. . Reported are: 


Month and Day: Date event occurred. 


Begin, Maximum and End: Begin, maximum and end times of 
the x-ray burst. The begin time is often well defined 
due to the impulsive beginning of most flare associated 
X-ray events. In the past, x-ray events were ended 
when the flux levels decayed to a point halfway between 
the maximum level and the pre-flare background level. 
This definition occasionally resulted in unterninated 
events. For operational purposes the end time of an x- 
ray burst is now determined by both the peak flux and the 
decay rate. An X-class event is ended when flux levels 
decay to less than M5 levels and remain stable for a 
minimum of four minutes or when flux levels decay to less 
than the Ml level, regardless of stability. All bursts 
smaller than X-class are now ended when flux levels decay 
to the C5 level or when flux levels decay to a point 
halfway between the maximum level and the »re-flare 
background level and remain stable for a zinimum of four 
minutes. The typical x-ray signature dictates that the 
end time of a burst will be less well defined than the 
begin time. 


X-ray Class: As defined earlier in this descriptive text. 


OPT: Optical information (iÉ available) as follows: 
Type of Activity: FLA - Solar flare 
EPL - Eruptive prominence on the 
solar limb 
BSL - Bright surge on the limb 
587 =- Spray (including ejecta) 
LPS - Loops 
DSF - Disappearing solar filament 
UNK - Unknown activity 
Importance: S = Subflare 
area < 2.0 square degrees 
1 = Importance 


2.0 area 5.1 square degrees 


2 = Importance 


5.2 < area 12.4 square degrees 


U JAN ja A 


3 = Importance 
12.5 < area < 24.7 square degrees 


4 = Importance 4 
24.8 < area in square degrees 
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Brightness qualifiers F, N or B are generally appended to 
the importance character to indicate Faint, Normal or 

Е2 1 рин Ы сс угу, 

Location: Flare location in heliographic coordinates. 
Region: SESC-assigned region number. 


Radio: Flux unit peak value above pre-burst background 
of associated radio bursts at frequencies of 245, 
415, 606, 1415, 2695, 8800 and 13400 MHz. 
1 Flux Unit = 10° z watt m 2 Hz 1. 


Sweep: Sweep frequency radio event associated with the energetic 
event, as UE (generally only II, IV and CIM are reported: 
Type: I Noise storm | 
LI ow drifte burst 
1114-5505. diTf£ burst 
IV - Broad-band smooth continuum burst 
V - Brief continuum burst generally associated 
with Type LII bursts 
CIM - Long-lived continuum associated with storms 
of Type IV bursts 


SIM -. Storm.of Type ITIL bursts 
MWB - Short microwave continuum 
SDF - Decimeter slow drift burst 


Importance: Relative importance on a scale of 1 to 2. 
An importance 1 sweep event would be ccnsidered 
small or minor. An importance 3 event is a 
major or significant event. 


SID: Sudden Ionospheric Disturbance information, 
as follows: 


Type: SWF - Short-wave fadeout Х 
Frequency: Highest frequency SWF was detected. 
Importance: Subjective measurement of the relative 


importance of the event on a scale of 
1 to 3 (much the same as already described 
for sweep frequency radio events). 


Flare List 


A listing of solar flares (optical importance and/or x-ray class) 
reported during the previous seven days. The list has been edited 

on a daily basis by the duty solar forecaster to delete multiple 
reports of the same event. Criteria used in determining which 

report, or combination of reports, will be used include: Quality of 
the observation (observing conditions, was flare observed in its 
entirety?), type of report (real time or film review), and the 
consensus of reporting OP oe a ee: The List includes the following 


information: 
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Month and Day: Date event occurred. 


Begin, Maximum and End: Begin, maximum and end times 
erschesoperee u: are. If tche start of che flare was not 
observed, a "B" is appended to the begin time indicatinz 
the flare started "Before" that time. If the maximum was 
not observed, the space is left blank. If the end of tke 
flare vas not observed, an "A" is appended to the end 
time indicating the flare ended "After" the end time. 

Location: Flare location in heliographic coordinates. 

IMP: Importance of the optical flare, as defined earlier 
in this descriptive text. 

REG: SESC-assigned region number or left blank for unassociatec 
xeray events and for flares which occasionally occur in 
undefined regions. 

Class: Denotes x-ray class of the event according to 
the C, M and X classification described earlier in 
this descriptive text. 


Weekly Region Summary 


A summary of all regions visible on the solar disk during the 
previous seven days. Each date's assessment of a region is based 
upon observations made during the 24-hour UT day aná the location 
is corrected to 24002 on the date indicated. 


This summary information is compiled from the Joint AFGWC/SESC 

Solar Region and Activity Summary which is a composite picture of 
each region from all sunspot reports received at Boulder during the 
UT Day. The Weekly Region Summary contains the following information: 


MO and DA: Month and day 


Reslen Location: Location of the ‘midpoint of the latitudinal 
and longitudinal extremeties" of the white-light sunspot 
group in heliographic coordinates at 2400Z on the date 
reported. The position of the center of H-alpha plage 
associated with the region is given if sunspots were nct 
visible. 


Sunspot: Sunspot information is as follows: 
Area: Total area of the region's sunspots in millionths 
of the solar hemisphere. 
2: Modified Zurich (McIntosh) three letter (Zpc) 
classification as follows: 


Z = Modified Zurich Brunner classes, A through H 
p = Penumbra type of largest spot in group: 
x - no penumbra, r - rudimentary, s - small 
symmetric, a - small asymmetric, h - large 
symmetric, and k - large asymmetric 
c = Relative sunspot distribution or compactness of 
group: x - single spot, o - open, i - intermediate, 
and e - compact 


98 








BR. “Total number of individual sunspots in the group or 
region. 
MAG: Magnetic classification, as shown below: 
A = Alpha (single polarity spoc) 
B = Beta (bipolar spot configuration) 
G Gamma (atypical mixture of polarities) 
BG = Beta-Gamma (mixture of polarities in a dominantly 
bipolar configuration) 


BD = Beta with Delta configuration 
BGD = Beta Gamma with a Delta configuration 
D = Delta (opposite polarity umbrae within single 


Penumbra) 
H-Alpha Plage: Indicates no sunspots visible but 
region still visible in H-alpha spectro- 
heliograms. 


Calcium: Calcium spectroheliogram data from Manila, Philippines, 
as follows: 
Area: Total region area of calcium plage in millionths 
of the solar hemisphere. 
BRT: Calcium plage brightness on a subjective scale of 
107 1ov to 5.0.Chien):. 
Blank space indicates no calcium plage report available 
тәсиси region on that date. 


Flares: Tabulation of the number of C, M and X-Class x-ray 
bursts, number of subflares and number of importance 1, 2 
and 3 flares observed on that date. 


At the bottom, the average Carrington longitude of the region and the 
x-ray burst and flare totals are listed. Finally, a comment on ‘the 
fate of the region is included. This indicates whether the region 
crossed over the solar limb, died on the solar disk, or was still 

an active region on the disk at the end of the reporting period. 


Boulder Geomagnetic Substorm Log 


This is a tabulation of substorm occurrences as observed by ground 
magnetometers in the North American International Magnetospheric 
Studies (IMS) network. This network is comprised of twenty-six 
Stations; extending from north of the auroral oval to the equatorial 
zone. A substorm is a localized geomagnetic disturbance which 
usually occurs near local midnight and is restricted in longitude. 
However, the current systems developed during a substorm may affect 
ground magnetometers in the entire nighttime sector. Additionally, 
individual substorms may occur at local times, away from midnight, 
and may be as large as 24 hours (global) in longitudinal extent. 

By noting the time, location and scale of a substora; one may 
estimate the effect of a substorm at a specific location. Among 
the many substorm effects are ionospheric effects (which influence 
radio communications) and telluric effects (which may disturb long 
distance electric power and communications systems). The familiar 
aurora is a visible manifestation of the geomagnetic substorm. 
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Currently, the Log provides the date, onset time (UT) and direction 
Eon Boulder) of each substorm. The direction is listed as "east", 
"west" or "centered" (at Boulder's geomagnetic longitude). The 
comment section further describes the gecmagnetic field for a 
Particular day. Occurrences such as sudden commencements (sc), 
sudden impulses, magnetic storm periods, multiple substorms and 
substorm (relative) strengths are included. 


Recent Monthly Solar Indices (Preliminary) 


This is a listing of monthly preliminary solar indices for the 
past year (updated within one week after the end of each month). 
The daily sunspot numbers appearing on the first page of the 
Report, which are a consensus of observations from es many as half 
a dozen observations, are simply averaged and is reported as "SESC 
Observed". 


irreal Zurich sunspot numbers are reported as "Zurich" as they 
become available. The ratio of Zurich to SESC sunspot numbers are 
listed under "Ratio Zurich/SESC'"; these can be used to determine 
2557 :: k (correction) factor. In order to provide the best 
early estimate of the Zurich sunspot number, e.g. near the end oz 
the month when the number is not yet available, an SESC corrected 
number will appear in parentheses in the "Zurich" column. 


The monthly mean value of 10.7 centimeter radio flux (1 flux unit 
= 10 22 watt m 2 Hz l) is listed for comparison with monthly meen 
sunspot numbers for comparable periods. These radio flux values 
have been adjusted to 1 A.U., except where values appear in paren- 
theses. Parenthetical values are the average of the daily values 
transmitted to SESC in real time. Adjusted values are those which 
appear in Solar Geophysical Data (SCD). 


Special pages 


To some extent the Preliminary Report is regarded as an experimental 
publication by the NOAA/SEL Space Environment Services Center. “he 
Center publishes new data, often on a trial basis; within the 
following limitations: 


JA Mere be reason to believe that new data will 
be of interest to a broad range of Report subscribers. 
2 mels ion óf new data Must not have a substantial 


impact on the costs of producing the Report (as 
would, for example, publication of photographs). 


Some items added to the Report are of a "one-time-only" nature. 
Among these are special drawings, descriptions or reports of unusual 
solar-geophysical features or events. Other data forms, such as 

the tabulation of monthly mean sunspot numbers, are candidates for 
permanent inclusion in the Report. Because a full discussion of 
special items is not possible in this descriptive text; explanations 
will accompany new data forms. Subscribers are encouraged to 
comment on the utility of various new data provided in the Report. 
Suggestions regarding new items to be placed in the Report will be 
carefully considered. 
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SMS/GOES X-ray and Particle Plots 


The Synchronous Metevrological Satellites (SMS-1 launched in June 
1974 and SMS-2 launched in February 1975) and the Geostationary 
Operational Environmental Satellites (GOES-1 launched in October 
2757 250ES=-2 launched in June 1977 and GOES-3 launched in June 
1978) carry 0.5 - 4A and 1 - 8À solar x-ray detectors and seven 
Proton sensors covering the range 0.8 to 500 MeV. These sensors, 
as well as a number of other sensors, are currently monitored 
continuously (except during eclipse periods) in real time bv the 
NOAA/SEL Space Environment Services Center in Boulder.  GOES-2 and 
GOES-3, in geostationary orbis at 110° West and 135° West, respec- 
tively, are now in use; and SMS-1l and -2 and GOES-1 are in standby 
mode. 


Preliminary graphs of SMS/GOES x-ray and vroton data are printed in 
the Report as solar activity warrants, with noteworthy x-ray events 
being generally more frequent than solar protons events. 


The top line of the display identifies the satellite, the sensors, 
and the traces, as follows: 


Satellite: Either GOES-1 or GOES-2 
Sensors: Either XRA (x-ray) or EPS (energetic particles) 
Traces: XRA (line) 1 - 8 

A n (points) 0.5 - 4À 


EPS (line) Flux » 10 MeV 


The second line indicates the date and time of the start of the 
plot. The number appearing in parentheses, either (1) or (5) 
indicates whether the plots are based on l-minute or 5-minute 
averages of the data. 


The ordinate units of the XRA plots are in watt п72 and Che согге- 
sponding C, M and X x-ray classification threshold is labelled on 
pues margin of the plot. The ordinate units of the EPS plots 
2501011217 0” sec İ ster İ. The abscissa is a time scale in minutes 
or hours, indicating the time elapsed since the beginning of the 
plot at the specified start time. 


Bad data have been deleted vherever possible. 


Last page 


The last page contains the current preliminary H-alpha synoptic 
chart. These maps depict the large-scale distribution of solar 
magnetic fields as inferred from the positions of filaments (cross- 
hatched), plage corridors (solid lines through stippled areas), 
filament channels (solid lines outside active regions) and sunspots 
(large dots, drawn very schematically). The stippled areas represent 
the H-alpha plages. Dashed lines are interpolations end estimates 
required to obtain consistency with polarities and patterns observed 
on previous solar rotations. The numbers identify active regions. 
The "plus" and "minus" signs indicate the polarity of the solar 
Magnetic field. The heliographic latitude is given along the left- 
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and right-hand side of the map, the Carrington longitude is printed 
along the bottom, and the date of the central meridian passage of a 
pavwen longitude appears at the top. Some portions of the maps may 
be inaccurate; however, every effort is made to correct and refine 
the maps as more detailed study and comparison of time sequences of 
H-alpha photographs and solar magnetoerams permit. Edited charts 
for complete Carrington. rotations appear later in Solar Geophysical 
Data (SGD). 


Notes: 


[x All times used within the Preliminary Report are Universal 
Times or Greenwich Mean Times (as indicated by UT, GMT or 
2: 


2. The highliyhts on the first page and all daca listings 
cover the seven-day period Monday through Sunday. The 
27-day forecast period begins on the following Wednesday. 


ə. All data listings are prepared by computer. Occasionally, 
inaccurate data will be reported and included. SESC 
appreciates notification of any guestionable data. 


4. The following are observatories whose reports are included 
in the SESC data base and utilized for preparing the 
Preliminary Report and Forecast of Solar-Geophysical 
Activity: 
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STATION 

Anchorage, Alaska, USA 
Athens, Greece 

Boulder, Colorado, USA 
College, Alaska, USA 
Culgoora, Australia 

Ft. Yukon, Alaska, USA 
Fredericksburg, Virginia, USA 
GOES-2 (USA satellite) 


GOES-3 (USA satellite 


Hiraiso, Japan 


Holloman AFB, New Mexico, USA 
i kutsk, USSR 

Kitt Peak, Arizona, USA 
Learmonth, Australia 

1 Posta, California, USA 

Manila, Philippines 

McMath Obs., Pontiac, Michigan, USA 
Moscow, USSR 

MSFC (Huntsville), Alabama, USA 
Ottawa, 
Palehua, 


Poitiers, 


Canada 

Hawaii, USA 

France 

Ramey, Puerto Rico 

Sacramento Peak, New Mexico, USA 
Sagamore Hill, Massachusetts, USA 
Sheep Mountain, Alaska, USA 
Sydney, Australia 

Thule, 
Tokyo, Japan 


Greenland 


Vostok, Antarctica 
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WMO ID 


702739 
16716 


— gem =~ gem e 


emp em = wn, em 


URSI ID 


25601 
32401 
20401 
25602 
55303 


EE 


N61 
N38 
N40 
N64 
S30 
N67 
N38 


geostationary 


geostationary 


N36 
N33 
N52 
N32 
s 
N34 
N14 
N42 
N55 
N35 
N45 
N21 
N46 
N18 
N32 
N42 
N62 
S34 
N76 
N36 
S78 


W159 
E-23 
W105 
W147 
E149 
W145 
W 77 
W109 


UN 


El41 
W106 
E104 
Wlli 
Ell 
555: 
E 
W 83 
ET 
W 87 
W 76 
W158 
EMEN 
W 67 
W105 
W 70 
W148 
Elo: 
W 68 
E140 
Е105 





SUMMARY OF SERVICES FROM THE NCAA SPACE ENVIRONMENT SERVICES CENTER 


* 


Uierweeklyo publication Preliminary Report and Forecast of Solar 
and Geophysical Activity." A summary of activity indices end 


event data for the preceding week and a prediction for the next 


27 days. 


A tape recorded message of solar and geophysical activity and 
indices for the most recent 24 hours and for the next 24 hours. 
Teleprenemnenbers rare: Commercial: (303) 499-8129 or FTS: 
32323235, 


A 40-second message with similar information on WWV at 18 
minutes past each nour. 


EE EE end Geophysical Actisiswe issued each 

ay at 2200UT. Contains activity summary, solLar and geomagnetic 
activity indices tor previous 24 hours and predictions for the 
next 72 hours. Available by teletype with recipient expected 

to pay cost of transmission. 


HF Radio Propagation Report issued each day at 0600UT (with 


s eendary wer sazes ae UUOUUT, 1200UT, 1S00L%). Distribution is 
also by teletype with recipient expected to pay costs. 


Direct inquiries to the duty forecaster (0700 to midnight, seven 
ayS per weex, mountain Time). Call (305) 299-1000 extension 
or ror TeS We do not accept collect calls). 


iwestenotltrieation by phone or Teletype. Persons or organizaticrs 
e Mediate notification of the precictionsor occurrence 
of various solar geophysical phenomena are contacted when their 
alert thresholds are met. 


Direct access to the SELDADS (Space Environment Laboratory Data 
Acquisition and Display System). Using standard computer termines, 
customers may access the SELDADS and obtain printouts of solar, 
interplanetary and geomagnetic variations data. 


Descriptive leaflets of each of these are available. For more 
information, please write to: 


Space Environment Services Center 
NOAA R43 

325 Broadway 

Boulder. 80303, USA 


The Space Environment Services Center is a real time service 
organization for serving those people needing real time or 
near-real-time operational support. If you need archival data 
for research studies, .you should write to: 


National Geophysical and Solar-Terrestrial Data Center 
Environment Data Services 

D6 NOAA 

325 Broadway 

Boulder, Colerado 80303,7 USA 


M. 


104 





i 
d 
| 
i 
| 





APPENDIX D 


IRANSFER FUNCTION USAGE 


The spectrum analyser when in the dBV display mode 
reads: 
E A volts 
N (dBV) 22010g Lone 


0 dB referenced to l nT and not l volt is desired. Thus: 
A volts, (X zr (E volts we: volts ) 
IO Lt X volts/nT X volts/nT” “1 volt nT 


where B (the field measured by the system) is: 


( 


A volts 


B= (z am) 


so referenced to 1 nT we have: 





N(dB )22010g + 20log X 


B 
iar 
X is the system transfer function given in units of 
volts/nT. Since spectrum analysers use various averaging 
vindovs and have various numbers of bins the folloving 


correction must be applied: 


B2=% (BW) , 


and: 


(Bw) = (£requency range, window 


number of bins’ ‘correction 


Now we have: 
n(a8)=20108[01/2 (8w)1/23+20108x 


and finally after rearranging: 
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10logö - N(dB) - 20logX - 10log(BW) , 


where l0logó is the corrected reading or the power spectral 


density. 9 has the units of nT^/Hz. 
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